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FOREWORD 


This  report  was  prepared  in  the  Mathematics  Research  Branch,  Flight 
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ment Command  under  RDO  No.  460-50-25,  "Mathematical  Investigation  of  Pilot  Escape 
from  High  Speed  Aircraft."  Harry  Ferguson  acted  as  project  engineer. 
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ABSTRACT 


It  is  shown  that  the  maximum  acceleration  that  a person  may  encounxer 
at  a given  time  after  separation  of  his  escape  unit  from  a damaged,  unrecover- 
able aircraft  is  dependent  only  on  the  velocity  of  the  aircraft.  Acceleration 
histories  determined  from  flight  tests  indicate  that  these  theoretical  maximums 
are  reasonable  and  practical  at  least  in  the  subsonic  case. 

Lt.  Colonel  J.  P.  Stapp  of  the  Aero  Medical  Laboratory  has  reported 
(see  paragraph  3)  that  high  rate  of  acceleration  onset  during  a fixed  time 
produces  injuries  not  found  at  lower  onset  rates.  On  the  hypothesis  that 
decay  rates  may  be  of  comparable  interest  it  is  shown  in  this  report  that 
decay  acceleration  rates  are  functions  of  the  initial  acceleration  and  speed 
of  the  escape  unit. 

The  accelerations  parallel  to  and  normal  to  the  spine  of  a person  in  a 
tumbling  escape  unit  are  calculated  for  an  assumed  hypothetical  case.  The 
effects  of  these  alternately  increasing  and  decreasing  accelerations  on  the 
person  may  be  described  as  a shaking  phenomena  which  may  be  beneficial  to 
the  escaning  person. 

The  security  classification  of  the  title  of  this  report  is  UNCLASSIFIED. 


PUBLICATION  REVIEW 

This  report  has  been  reviewed  and  is  approved. 
FOR  THE  COMMANDING  GENERAL: 


Colonel,  USAF 

Chief,  Flight  Research  Laboratory 
Directorate  of  Research 
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LIST  OF  SYMBOLS 

= time  (seconds)  measured  from  the  instant  of  separation  of  escape  unit 
= distance  (feet)  travel  of  the  escape  unit 
= speed  and  acceleration,  respectively  after  t seconds 
= speed  and  acceleration  respectively  at  separation  time  t = o 
= effective  drag  area  of  the  escape  unit 

= drag  coefficient  of  escape  unit  (considered  constant  in  this  study) 

= weight  (pounds)  of  escape  unit 
= altitude  (feet)  above  sea  level 

o 

a density  (slugs/ft  ) of  atmosphere 
■M  = jolt  measure 


WADCTR  52-278 


v 


INTRODUCTION 


The  nature  of  the  accelerations  and  jolts  that  a pilot  may  expect  to 
encounter  during  escape  from  a damaged,  unrecoverable,  high  speed  aircraft 
would  appear  to  depend  on  many  conditions  such  as  vehicle  of  escape,  speed 
of  the  aircraft,  altitude  of  the  aircraft,  motion  during  escape,  mass  of 
escape  unit,  and  the  method  of  holding  the  pilot  in  the  escape  unit.  These 
conditions  and  the  physical  and  psychological  tolerances  of  humans  to  accel- 
erations determine  design  criteria  for  escape  mechanisms  and  associated 
aircraft . 

Although  all  the  factors  mentioned  above  are  not  treated  in  this  report, 
the  major  part  of  them  are  considered  and  the  relative  importance  of  these 
several  parameter  conditions  is  indicated.  This  study  does  not  consider 
the  process  or  the  effect  of  the  build  up  to  the  peak  acceleration. 
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Since  the  accelerations  of  primary  interest  will  he  of  short  duration 
(less  than  five  seconds)  and  of  orders  of  magnitude  above  2g,  the  gravitational 
force  component  which  in  general  will  be  essentially  normal  to  the  path  has  been 
ignored  in  this  study.  Also  no  attempt  is  made  to  take  into  account  the  effect 
of  pitching,  pitch  attitude  and  acceleration  on  the  drag  coefficient.  It  has 
been  anticipated  that  compressibility  will  not  cause  abrupt  changes  in  drag  at 
initiating  speeds  below  the  transonic  and  during  the  escape  period  of  interest. 
For  initiating  speeds  at  or  above  transonic  speed  a discontinuous  or  abrupt 
Change  in  drag  may  be  expected  and  therefore  the  results  are  subject  to 
correction  for  this  effect.  A subsequent  study  will  treat  this  matter.  Under 
the  conditions  specified  the  equation  of  motion  becomes. 

(1)  x = -K  x2 


in  which  K is  a positive  constant.  This  equation  has  solutions  defining  the 
time  histories  of  speed  and  acceleration  of  the  escape  unit.  These  are 


and 


(2)  _i_ 

xn 


= 2 

1 - *£  t 


(x  is  negative) 


(1 


7 


x 

o 


It  is  of  interest  to  note  how  these  histories  depend  on  the  initial  speeds  and 
accelerations.  The  decay  rates  of  the  speed  and  acceleration  depend  only  on 
the  ratio  xQ  and  the  latter  is  more  rapid  than  the  former. 


2.  Estimated  Upper  Bounds. 

As  has  been  indicated  in  the  introductory  remarks  a multiplicity  of  escape 
conditions  may  generally  be  anticipated  and  it  is  desirable  to  develop  overall 
information  which  depends  on  as  few  as  possible  of  these  parameters.  This  has 
been  accomplished  in  this  instance  by  calculating  the  maximum  acceleration  that 
may  be  expected  at  some  given  time  t after  initiation  of  escape. 

For  a given  t and  escape  speed  x0  the  maximum  acceleration  that  may  be 
expected  is  found  by  setting  the  derivative  of  x with  respect  to  x0  equal  to 
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zero  and  solving  tor  x*  in  terms  of  t and  x . 

o o 

("*.  - 0*1  - IxZtf 

/Kc  I'***  tCcX) 

{sX0  - x7CQtf 

and  thus  the  *x*  which  provides  this  is  evidently 

(5)  = 

t 

and  from  (3)  the  maximum  acceleration  takes  the  form 

(6)  maxjx/^x^ 

4t 

Clearly  this  result  is  independent  of  weight  of  escape  unit,  drag  coefficient, 
effective^drag  area,  and  altitude  of  escape.  From  (5)  it  is  also  clear  that 
the  max  I x I is  always  attained  at  some  time  during  escape  for  each  initial 
escape  speed. 

Equation  (6)  thus  provides  upper  bounds  for  acceleration  histories  during 
escape  which  depend  only  on  the  initial  velocity  of  escape.  Of  course  these 
upper  bounds  are  not  practical  for  values  of  t close  to  zero  and  an  alternate 
method  is  used  for  these  small  time  values.  Figures  1,  2,  and  3 show  graphs 
of  (6)  for  the  time  interval  0.5  to  4-0  seconds  and  for  initial  speeds  at 
200  MPH  intervals,  from  400  MPH  to  1800  MPH. 

Figure  1 shows  a pair  of  acceleration  histories  determined  from  flight 
tests.  These  curves  when  compared  with  the  theoretical  curves  indicate  that 
the  latter  are  reasonable  and  practical  upper,  bounds  in  the  subsonic  case 
of  acceleration  histories  over  a wide  range  of  escape  conditions  (for  more 
detail  see  paragraph  4 below).  It  should  be  noted  however  that  these  curves 
include  speeds  which  are  subsonic,  transonic,  and  supersonic  and  tint  the 
equation  of  motion  (1)  is  not  completely  applicable  for  motions  initiated 
at  supersonic  and  transonic  speeds.  Attempts  to  obtain  more  quantitative 
information  in  these  cases  are  in  process. 

The  following  method  was  used  to  establish  upper  acceleration  bounds  in 
the  time  interval  0.0  to  0.5.  At  t = o we  have 


We  have 


• # 


JLJL 
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(7)  x'0  = - gSCpP^o2 

2W 

Examination  of  experimental  results  indicate  that  l/30  may  be  used  for  an  upper 
bound  for  S CD  and  hence 

W 

(8)  i*0  _i_  x0  2g 

60 


The  upper  bound  at  t=o  is  estimated  by  (8)  for  a specified  altitude  and  was 
joined  by  a straight  line  segment  to  the  upper  bound  at  t=0.5  estimated  by  (6). 
This  line  segment  determines  conservatively  the  upper  bound  for  the  time  interval 

0. 0  sL  t jsi  0.5  and  for  the  specified  altitude  as  a minimum  altitude.  Since 
this  result  depends  on  the  density  altitude  the  results  as  presented  in  figures 

1,  2,  and  3 correspond  in  this  time  interval  0.0  to  0.5  seconds  to  the  minimum 
altitudes  there  indicated. 


When  using  the  graphs  in  figures  1-3  it  is  important  to  bear  in  mind  that 
the  acceleration  man  can  tolerate  is  dependent  on  the  time  interval  and  the 
specific  attitude  of  the  man  during  the  acceleration.  Man  can  tolerate  higher 
acceleration  for  shorter  time  intervals,  and  can  endure  greater  accelerations 
in  direct ions. normal  to  the  spine  than  in  directions  parallel  to  the  spine.  In 
the  direction  parallel  to  the  spine,  man  can  stand  greater  acceleration  in  the 
direction  toward  the  head  (blood  flows  from  the  head  causing  black  out)  than  he 
can  stand  in  the  direction  away  from  the  head  (blood  flows  to  the  head  and  may 
cause  red  out). 

3.  The  Jolt. 


It  has  been  shown  that  the  rate  of  onset  of  acceleration  can  produce 
serious  injury  to  the  submitting  subject  as  reported  in  AF  Technical  Report 
No.  5915,  "Human  Exposures  to  Linear  Deceleration,  Part  2:  The  Forward- 
Facing  Position  and  the  Development  of  a Crash  Harness"  by  Major  John  Paul 
Stapp.  This  effect  is  defined  as  a jolt  and  the  third  derivative  "x  has  been 
adopted  as  its  measure.  Using  the  same  measure  for  deceleration  decay  effect 
and  again  assuming  that  equation  (1)  is  applicable,  that  is,  that  the  drag 
force  is  proportional  to  the  square  of  the  speed,  one  finds  readily  that 


Since  x0 


(9)  x*  = 2*o2*o2 

(xQ  - x’ot)3 

is  negative  it  is  obvious  that  x*  is  a maximum  when  t=o  and  hence 


(10) 


Figure  (4)  shows  graphs  of  (10).  The  graphs  show  the  jolt  as  a function  of  the 
initial  acceleration,  for  initial  speeds  at  200  MPH  intervals  for  a range  from 
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400  MPH  to  1800  MPH.  The  graphs  for  initial  velocities  of  400,  600  and  600  MPH 
were  cut  off  at  values  of  x’Q  which  correspond  to  the  upper  bounds  of  x estimated 
at  t=0.  The  other  curves  would  have  correspondingly  been  cut  off  if  the  x0 
scale  had  extended  sufficiently  far  beyond  76g. 

It  should  be  noted  that,  contrary  to  the  upper  bound  results,  jolt  results 
are  probably  valid  at  both  the  transonic  and  supersonic  speeds.  On  the  other 
hand  they  measure  only  the  jolt  due  to  acceleration  decay  and  give  no  information 
about  onset  effects. 

• • 

It  is  interesting  to  note  that  for  a given  initial  acceleration  xc  condition 
the  jolt  decreases  with  increase  in  initial  speed  and  that  for  a given  initial 
speed  the  jolt  increases  with  increase  in  initial  acceleration.  Thus  an  aero- 
dynamic clean  up  of  the  escape  unit  vail  reduce  jolt  while  a reduction  in 
weight  will  increase  jolt. 

4.  Experimental  Data  Shows  Maximum  Decelerations  and  Jolts  are  Nearly 

Attained. 


An  investigation  was  made  to  determine  the  accelerations  encountered  in  tests 
of  seat  ejections  carried  on  at  Muroc  and  Hamilton  Air  Force  Bases.  These  tests 
are  reported  in  Memorandum  Report,  MR#MCREXA7-45341-4-*l,  Pilot  Ejection  Flight 
Tests  Conducted  with  a TF-80C  Airplane  at  Muroc  and  Hamilton  Air  Force  Bases. 
Graphs  of  the  paths  of  the  ejections  are  given  in  the  report.  Using  the  graphs 
of  Test  No.  22,  Page  27,  live  test  (Captain  V.  Mazza),  and  Test  No.  36,  Page  41, 
dummy  test,  horizontal  and  vertical  accelerations  were  determined  mathematically. 
The  graphs  of  these  are  plotted  on  figure  1,  and  the  maximum  jolts  encountered 
are  marked  on  figure  4.  It  is  observed  that  the  experimental  values  of  accelera- 
tion and  jolt  attain  values  very  near  the  theoretical  maximums. 


5.  Explanation  of  the  Graphs  of  x as  a Function  of  *x*„. 

« • ••  * 2 

Graphs  (figures  5-7)  have  been  drawn  using  x = x°x°  for  t = 

. » #* 77 

(Xo  - Xot) 

1,  2,  3,  and  xQ  = 400,  800  and  1200  MPH  and  letting  x0  vary.  It  is  of 

interest  to  note  that  for  a given  time  t,  a given  velocity  of  the  aircraft,  and 
for  a given  acceleration  of  the  man  there  are  two  initial  accelerations.  For 
example,  figure  6,  if  x0  = 800  MPH  there  are  two  initial  accelerations  -8g 
and  -38g  which  give  an  acceleration  of  -4g  when  t = 2.  The  relationship 
between  values  of  time,  velocity  of  the  aircraft,  and  the  two  initial  accelera- 
tions x0  and  x00  is  obtained  by  solving 


**  * 2 
xoxo 


U0-x0T)2 


* * *2 
x00x0 


(x o-Xoo^) 2 
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400  MPH  to  1800  MPH.  The  graphs  for  initial  velocities  of  400,  600  and  800  MPH 
were  cut  off  at  values  of  x’0  which  correspond  to  the  upper  bounds  of  x estimated 
at  t=0.  The  other  curves  would  have  correspondingly  been  cut  off  if  the  x0 
scale  had  extended  sufficiently  far  beyond  76g. 

It  should  be  noted  that,  contrary  to  the  upper  bound  results,  jolt  results 
are  probably  valid  at  both  the  transonic  and  supersonic  speeds.  On  the  other 
hand  they  measure  only  the  jolt  due  to  acceleration  decay  and  give  no  information 
about  onset  effects. 

It  is  interesting  to  note  that  for  a given  initial  acceleration  xQ  condition 
the  jolt  decreases  with  increase  in  initial  speed  and  that  for  a given  initial 
speed  the  jolt  increases  with  increase  in  initial  acceleration.  Thus  an  aero- 
dynamic clean  up  of  the  escape  unit  will  reduce  jolt  while  a reduction  in 
weight  will  increase  jolt. 

4.  Experimental  Data  Shows  Maximum  Decelerations  and  Jolts  are  Nearly 
Attained. 

An  investigation  was  made  to  determine  the  accelerations  encountered  in  tests 
of  seat  ejections  carried  on  at  Muroc  and  Hamilton  Air  Force  Bases.  These  tests 
are  reported  in  Memorandum  Report,  MR#MCREXA7-45341-4-l,  Pilot  Ejection  Flight 
Tests  Conducted  with  a TF-80C  Airplane  at  Muroc  and  Hamilton  Air  Force  Bases. 
Graphs  of  the  paths  of  the  ejections  are  given  in  the  report.  Using  the  graphs 
of  Test  No.  22,  Page  27,  live  test  (Captain  V.  Mazza),  and  Test  No.  36,  Page  41, 
dummy  test,  horizontal  and  vertical  accelerations  were  determined  mathematically. 
The  graphs  of  these  are  plotted  on  figure  1,  and  the  maximum  jolts  encountered 
are  marked  on  figure  4.  It  is  observed  that  the  experimental  values  of  accelera- 
tion and  jolt  attain  values  very  near  the  theoretical  maximums. 

5.  Explanation  of  the  Graphs  of  x as  a Function  of  x'^ . 

« ® ••  • 2 

Graphs  (figures  5-7)  have  been  drawn  using  x = x°x°  for  t = 

(x0  - x0t) 

1,  2,  3,  and  xQ  = 400,  800  and  1200  MPH  and  letting  x0  vary.  It  is  of 

interest  to  note  that  for  a given  time  t,  a given  velocity  of  the  aircraft,  and 
for  a given  acceleratfon  of  the  man  there  are  two  initial  accelerations.  For 
example,  figure  6,  if  x0  = 800  MPH  there  are  two  initial  accelerations  -8g 
and  -38g  which  give  an  acceleration  of  -4g  when  t = 2.  The  relationship 
between  values  of  time,  velocity  of  the  aircraft,  and  the  two  initial  accelera- 
tions xc  and  x00  is  obtained  by  solving 

..  . 2 
*0*0 

(*0-x0T)2 
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for  t.  One  obtains  t 
are  equal: 


/ 


TT-n 

xoxoo 


where  t is  the  time  when  both  accelerations 


From  the  above  information  and  graphs  one  concludes  that  for  two  ejections 
A and  B at  same  initial  velocity  and  initial  accelerations  xQ,  x respectively, 

where  xQ  is  greater  than  xOD,  that  at  some  instant  t the  accelerations  are  equal 
and  prior  to  t the  acceleration  of  A is  greater  than  that  of  B and  after  t the 
acceleration  of  B is  greater  than  that  of  A. 


6.  Acceleration  .Estimates  for  a . Tumbling  Escape  Unit. 


Some  preliminary  acceleration  calculations  on  a tumbling  escape  unit  have 
been  made,  based  on  simplifying  assumptions  concerning  the  nature  of  the  tumbling 
and  the  drag  coefficient  during  tumbling.  Figure  I below  defines  the  basic 
geometry  considered. 


Fig.  I 

The  vector  labeled  x defines  the  displacement  of  the  c.g.  of  the  escape 
unit  in  the  direction  of  the  airstream.  The  angle  9 is  the  angle  between  the 
person^  spine  line  and  the  direction  of  motion,  and  w=9  is  the  angular  velocity 
of  the  persons  spine  line.  T is  the  distance  from  the  persons  c.g.  to  the  body 
point  at  which  accelerations  are  being  considered.  F^  and  F^  are  respectively 
the  acceleration  components  normal  and  parallel  to  the  persons  spine  line  of 
the  acceleration  at  some  given  point  in  the  body. 

The  tumbling  history  of  the  escape  unit  was  arbitrarily  defined  to  be 
w = wQ  (1  - e-’^)  in  which  wD  is  the  terminal  angular  velocity.  The  graph 
of  the  function  is  below. 
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Accelerations  fn*.  ft  each,  consist  of  the  sum  of  three  accelerations 
due  to  linear  accelerations  of  the  person,  rotation  of  the  person,  and  the 
force##of  gravity  on  the  person.  The  linear  acceleration  of  the  person  is 


x = 


xoxo 


so  its  contribution  to  F„  is  x cos  6 and  to  Fm  is  - x sin  9 

N JL 


(io-x'ot)* 

It  is  assumed  that  accelerations  in  direction  of  stomach  to  head,  and  in  direction 

-.6t 


The  rotational  motion  contributes  to  FN  an 


of  back 'to  chest  are  positive. 

acceleration  -w  y = -w02(i-e  and  to  F^  an  acceleration  of  wr  = .6w07'e 


Force  of  gravity  gives  an  acceleration  of  g sin  9 to  F„  and  g cos  9 to  F, 
Therefore  w 

Fjj  = x cos  0 -w2f'  +g  sin  9 

F =-x  sin  9 +wr+g  cos  9 
T 


T* 


Graphs  have  been  drawn  utilizing  these  formulas  (figures  8-25).  These  graphs 
were  drawn  using  the  following  values  of  the  parameters: 


x0  = 587(400  M.P.H.),  x0  = -4f0  . , T = 1,  2 and, 

wc  = J^-ff(100  r.p.m. ),  5_2l50  r.pm.),  j5Jt  (25  r.p.m.) 

J 3 6 

Values  of  r = 1,  2 are  selected  as  representative  of  distance  from  center  of 
rotation  to  the  heart,  and  brain  respectively  of  the  person.  Graphs  are  drawn 
of  F^figures  8,  11,  14,  17,  20,  23)  and  F^,  (figures  9,  12,  15,  18,  21,  24)  as 

functions  of  t.  To  construct  these  graphs  it  is  necessary  to  express  9 as  a 
function  of  t . 


& = w 
dt 

-,6t 

d9  = w0  (l-e  ) dt 


Integrate  and  use  initial  condition  that  9 = 9C  when  t=0  and  obtain 

® - «o  * wo  [t  - 1 U-e‘:6t)] 
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Graphs  shown  in  figures  10,  13,  16,  19,  22  and  25  show  the  actual  vector 
acceleration  that  the  tumbling  escaping  person  would  encounter.  The  effect 
might  be  described  briefly  as  a shaking  phenomenon.  The  damage  that  such  a 
phenomenon  might  cause  would  be  related  to  the  natural  frequencies  of  the  body 
part  concerned,  the  duration  of  the  shake,  and  the  frequency  and  amplitude  of 
the  shake. 
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